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ABSTRACT 


The dynamic response of the oceanic mixed-layer to the thermodynamic forcing 
at the sea surface is analyzed in order to describe the pattern of the oceanic meridional 
overturning. The technique proposed in this study is based on residual-mean theory, 
which takes into account the transport of buoyancy and tracers by transient eddies. From 
the observed air-sea density flux and mixed-layer density distributions, we estimate the 
two components of the Meridional Overturning Circulation (MOC) corresponding to the 
adiabatic (along-isopycnal) advection and the diabatic (cross-isopycnal) flux. 
Calculations are performed for the global ocean and, additionally, for each oceanic basin. 
The derived method extends the Walin (1982) water mass transformation theory, and 
permits, for the first time, assessment of the strength of the MOC adiabatic component 
from the sea surface data. 

This study offers a statistical description of the atmospheric and oceanic 
databases and gives some suggestions for the choice of specific datasets. In particular, the 
two most reliable atmospheric climatology databases (ECMWF and NCEP/NCAR re¬ 
analyses) are compared, and the impact of their inaccuracies on the MOC calculations is 
evaluated. 

The results presented in this paper are consistent with the pattern of the 
thermohaline circulation estimated from in-situ measurements and models. They also 
support the previous estimates of diapycnal volume flux by the Walin-type calculations. 
Eurthermore, they suggest that the global thermohaline circulation is dominated by the 
adiabatic advection in the ocean interior. 

One of the goals of the data analysis and numerical modeling effort is to explain 
the role of the mesoscale variability in the dynamic coupling of the ocean and the 
atmosphere. Properties of the ocean thermal fronts and eddies determine the Undersea 
Warfare (USW) tactics in the areas of high mesoscale activity, and therefore our efforts 
to predict their distribution and strength are directly related to the Navy research interests. 
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I. 


INTRODUCTION 


The ocean plays a major role in the climate system since it is the main absorber of 
solar energy and therefore constitutes the direct source of most of the energy in the 
atmosphere. Thus, a complete understanding of the complex interactions and feedbacks 
between the ocean and the atmosphere is essential for forecasting any climate variation. 

The interconnections between the heat and moisture fluxes at the sea surface and 
the oceanic and atmospheric circulations determine the equilibrium state of the climate 
system. Of particular importance are the meridional transports of energy by winds and 
currents (wind- and density-driven). As a consequence, fluctuations in the density-driven 
oceanic Meridional Overturning Circulation (MOC), also referred to as thermohaline 
circulation, may lead to major climate variations. The MOC has recently become a 
subject of great interest in the context of the currently observed global warming in polar 
regions. 

One of the biggest uncertainties in the theory of the MOC is associated with the 
relative importance of the interior mixing and the adiabatic advection for its 
establishment and maintenance. The classical theories use purely diffusive models to 
depict the thermohaline circulation: Stommel (1961) emphasized the importance of the 
thermodynamic forcing on the individual density components, and Munk (1966) 
described the vertical mixing as the means for communicating the buoyancy anomaly 
from the surface into the ocean interior. Even more recent comprehensive discussions 
(Welander 1986, and Whitehead 1995) upheld such a description of the MOC. 

However, several new theoretical and observational findings cast some doubt on 
the ability of these purely diffusive models to explain the magnitude and pattern of the 
meridional overturning. Indeed, although numerical simulations of the Atlantic exhibit a 
dependence of the MOC on diapycnal diffusivity (Bryan 1987, among others), most 
diffusive models have to rely on very high values of vertical mixing (5-10xl0^^m^5 ') 
to reproduce oceanic overturning of realistic strength. Yet these high values are not 
supported by the microstructure and tracer dispersion measurements which instead 
suggest the thermocline diffusivity on the order of lxl0^^m^5 ' (Ledwell et al. 1993; 
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Toole et al. 1994). This has led to an alternative adiabatic view of the meridional 
overturning, which emphasizes the importance of the advection along isopycnal surfaces, 
induced by wind forcing and eddy transfer (Toggvieler and Samuels 1998; Gnanadesikan 
1999; Marshall and Radko 2003; Radko 2006). The relative importance of the interior 
mixing and adiabatic advection for the MOC establishment and maintenance has become 
the subject of a vigorous ongoing debate. 

The classical representation of the MOC, shown in Fig. 1 and Fig. 2, focuses on 
the abyssal circulation driven by a few deepwater sources where deep convection occurs, 
and a large-scale upwelling that compensates for the localized downward flux (Stommel 
1961). The modem theory of meridional overturning (Webb and Suginohara 2001; 
Bocaletti et al. 2005; among others) identifies and describes at least two distinct dynamic 
components of circulation, which are the shallow overturning cells in the main 
thermocline and the deep circulation in the abyssal ocean. In the abyssal regions, shielded 
from the direct influence of wind, small-scale mixing processes are necessary to re¬ 
supply the potential energy removed in the interior by the overturning and eddy¬ 
generating process (Wunsh and Ferrari 2004). On the other hand, the dynamics of the 
strongly stratified thermocline and the associated shallow overturning are controlled by 
the ventilation of water masses along the isopycnals that outcrop at the sea-surface 
(Luyten et al., 1983); therefore, the role of diapycnal mixing in the maintenance of the 
upper cell is questionable (Marshall et al. 2002). Indeed, unlike the abyssal regions, the 
thermocline is directly influenced by the surface winds, which provide the energy for the 
meridional overturning and make it possible to close the circulation even in the absence 
of vertical diffusion below the mixed-layer. 

Traditional views on the oceanic overturning tended to emphasize the dynamics 

of the deep overturning circulation (Roemmich and Wunsh 1985; Munk and Wunsh 

1998) and its major role in the meridional heat transport. However, the relative 

significance of the upper and deep circulation cells should be re-evaluated in the context 

of the recent suggestion (Boccaletti et al. 2005) that the meridional transport of heat is 

controlled by the processes operating in the upper ocean. In an attempt to meaningfully 

quantify partitioning of the oceanic heat flux between the deep and shallow overturning 

cells, these authors introduced the “heatfunction” which identifies the components of 
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circulation that effectively contribute to the total meridional heat flux. As shown in Fig. 
3, diagnostics of a numerical model (Boccaletti et al. 2005) have pointed out that the 
heatfunction maxima are associated with a surface intensified flow largely limited to the 
central thermocline. This suggests that the meridional heat transport is dominated by the 
contribution from the upper branch of the MOC, the dynamics of which are still poorly 
understood (Samelson 2004). 


heat loss HERE input 

from ocean to ocean 


tat t 1 M 



Figure 1. Schematic of the MOC (from http://www.liv.ac.uk/physocean/research.htm). 


Although the global characteristics of the oceanic overturning can be estimated 
from the air-sea fluxes (Trenberth and Caron 2001), the vertical pattern of the MOC is 
difficult to reconstruct from observations or models. Sparse oceanographic measurements 
allow only a crude partitioning of the MOC in terms of the volume fluxes in the shallow, 
intermediate, and deep waters (Talley 2003). All numerical modeling studies rely on the 
assumed parameterizations of diapycnal fluxes, and the magnitude and functional 
dependencies of these fluxes are still a greatest source of uncertainty. This difficulty, for 
instance, precludes a reliable quantitative estimate of the relative importance of the 
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adiabatic advection and cross-isopycnal diffusion for maintenance of the MOC. A 
successful approach in analyzing the meridional overturning was proposed by Walin 
(1982), who deduced the integrated rates of water mass transformation at isopycnal 
surfaces from the air-sea buoyancy fluxes and surface density data. Walin’s technique 
was extended and applied to the Atlantic (Speer and Tziperman 1992; Donners et al. 
2005) and the global ocean (Speer et al. 1995). The water mass transformation at a given 
buoyancy surface is directly linked to the cross-isopycnal mass flux (Tziperman, 1986), 
and this provides an estimate of the strength of the diapycnal (and therefore inherently 
diabatic) component of the oceanic circulation. 



A 


Figure 2. Classical representation of the global thermohaline circulation and its associated 


oceanic conveyor belt (from Aguado and Burt, Understanding Weather and 
Climate). 
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Assessing the adiabatic pole-to-pole component of the MOC is more challenging 
since the movement of fluid along the isopycnal surfaces does not significantly involve 
the interior water mass transformation. Another complication results from the presence 
of mesoscale eddies whose contribution to the overall overturning has to be taken into 
account. It is now widely accepted that the eddy transfer of buoyancy and tracers in the 
upper ocean is often comparable in magnitude to the advection by the time-mean flow 
(e.g. Radko and Marshall 2004; Henning and Vallis 2004). Since eddies generally tend 
to counteract the Eulerian transport, the residual flow (sum of the mean and eddy-induced 
circulations) is considerably different from the mean. Nevertheless, attempts have been 
made to diagnose the residual overturning circulation from the air-sea fluxes and surface 
buoyancy distribution (Speer et al. 2000; Karsten and Marshall 2002; Marshall and 
Radko 2003). These studies focused on the southern ocean, where predominantly zonal 
flow of the Antarctic Circumpolar Current justified the use of a simplified two- 
dimensional framework. A clear way of diagnosing the vertical structure of the global 
pole-to-pole MOC from the air-sea fluxes for complicated three-dimensional geometries 
of the isopycnals is still lacking. 

The purpose of this study is to combine Walin’s technique with the ideas of the 
residual circulation theory to quantify both the diabatic and the adiabatic components of 
the meridional overturning. More precisely, the along- and cross-isopycnal residual 
fluxes in the whole ocean and in each basin (Atlantic, Pacific, and Indian) will be 
estimated from the air-sea fluxes and mixed-layer buoyancy distribution. The technique 
presented here offers for the first time a quantitative assessment of the adiabatic 
component of the MOC which is based entirely on the sea surface data. 

The coordinate system introduced in this study is somewhat similar to the 
isentropic framework often used by meteorologists. Here, however, the isentrope field 
cannot represent the density distribution, since the seawater potential density also 
depends significantly on salinity. In addition, contrary to the isentropic circulation, the 
thermodynamic forcing in this MOC model is restricted to the sea-surface and the 
associated oceanic boundary layer. The circulation in the ocean interior is assumed to be 
a linear response to mixed-layer dynamics induced by the air-sea fluxes. 
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The formulation of the diagnostic technique is presented in section 11. Section III 
compares the different databases and suggests the choice of specific datasets. In section 
IV, the diagnostic model is applied to the global ocean and separately to the Atlantic, 
Pacific, and combined Indian-Pacific basins. Finally, the main results and the MOC 
characteristics are discussed in the conclusion (section V). 
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Figure 3. (a) Stream function related to the zonally-averaged flow inferred from the M.I.T. 

global circulation model. The distribution of meridional volume transport, given 
in Sv (10'’m^5 *), is defined by the stream function since seawater is nearly 
incompressible. Red shading corresponds to clockwise flow and blue to counter¬ 
clockwise flow. Contour interval is 5 Sv. (b) Zonal averages of the total heat 
transport of the global circulation model, (c) Meridional heat function in PW (10*^ 
W) inferred from the global circulation model. Red shading corresponds to 
northward heat transport and blue to southward heat transport. Contour interval is 
0.17 PW (from Boccaletti et al. 2005). 


6 






































II. FORMULATION 


The two surface properties used to reconstruct the large scale meridional 
overturning in the ocean are the mixed-layer density anomaly (cr^) and the air-sea 

density flux (B). The dimensionless mixed-layer density anomaly (cr^) is related to the 

mixed-layer density (in kg as follows: 


Pm 


= 1000 + o-„ 

m 


And the sea surface density flux (5) is conventionally defined as: 


^ aH . {E-P)S 
B = - + - 


C„ 


1-5 


( 1 ) 


where (a,/3) are the expansion/contraction coefficients of seawater, p^ is the standard 
seawater density, H is the heat flux from the atmosphere into the ocean, E and P are the 
evaporation and precipitation rates, 5 is the salinity, and is the specific heat capacity 

of water. The technique utilizes the observational data based on meteorological re¬ 
analysis (density flux at the sea surface) and Levitus database (seawater density) to 
reconstruct the pattern of the meridional overturning. 

For the computations in the Atlantic, Pacific, and combined Indian-Pacific 
oceans, we consider the density range where all isopycnals are zonally bounded by land, 
as sketched in Fig. 4. Confining the analysis to these water masses greatly simplifies the 
interpretation of the diagnostic results and accurately describes the upper ocean. For the 
global ocean, computations are extended to the maximum value of mixed-layer density 
anomaly. 


A. RESIDUAL FLUX AT THE BOTTOM OF THE MIXED-LAYER 

It has been shown (Andrews and McIntyre 1976) that the distribution of density 
and tracers in the eddying flows cannot be given by only the Eulerian mean circulation, 
but the distribution also involves the eddy-induced advection, a process similar to the 
Stokes drift. Therefore, the starting point is the time-mean density equation that takes 


7 



into account the eddy fluxes as follows: 


V ■VcT + V(yW) = — ( 2 ) 

dz 

where cr is the time-mean density anomaly, primes denote the perturbations from this 
mean due to transient eddies, and B represents the vertical density flux due to small-scale 
processes and air-sea fluxes. While we mostly associate the eddy fluxes with the 
mesoscale variability, it is convenient to include in it all transient processes, such as the 
effects of the seasonal cycle and variation in the mixed-layer depth. 

The eddy fluxes in Eq. (2) can be broken down into two distinct components 
(Radko and Marshall 2006): 

_ 

- The adiabatic advection by the eddy-induced velocity of the residual mean theory (V ), 
which is assumed to be non-divergent (V • V = 0); 


- The remaining diabatic component that can be written without loss of generality as the 
vertical convergence of the vertical flux 


'' as,,,/ 


dz 


As a result of this transformation, Eq. (2) becomes: 

t7 w 

V res-ya = — 

dz 


(3) 


where the residual velocity V res = V +V represents the advection of density and tracers 
by both mean field and adiabatic eddies, and B = B + B^^^^ includes the diabatic effects of 
the small-scale mixing and eddies. 

Eollowing Marshall and Radko (2003), we separately discuss the dynamics of the 
thin, vertically homogeneous mixed-layer ( -h^ < z < 0) and the stratified interior 

(z < -h^). In the mixed-layer the residual buoyancy equation reduces to: 


da da^ dB 

ftl ... HI 

^ ^ res ^ ^ 

ox oy oz 


(4) 


where a^ is the mixed-layer density anomaly. Hence, the integration of Eq. (4) over the 
depth of the mixed-layer results in: 
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V—+V — = B, + B,„-B 


Z=-h,„ 


(5) 


dx dy 

0 0 

with U = ^ i^resdz , y = j , Bq being the air-sea density flux, a eontribution 


-h^ 


-h,„ 


from the diabatie eddies in the mixed-layer, and B the vertieal density flux due to 

Z=-h,,. 


diabatie proeesses immediately below the mixed-layer. 

In the mixed-layer buoyaney equation (5), the total density flux on the right-hand 
side is elearly dominated by the direet air-sea foreing Bq . The estimates in Radko and 

Marshall (2006), in the eontext of the ACC, suggest that the sea surfaee flux B^ exeeeds 
^mie by an order of magnitude; the diabatie eddy fluxes play an even lower role on 
planetary seales (Radko 2006) that are the primary foeus of this study. In addition, both 
small-seale mixing and the diabatie eddy effeets are greatly redueed below the mixed- 
layer and therefore B.. » 5 . Henee, Eq.(5) may be approximated by: 




dx 


dy 


( 6 ) 


The integration of the residual eontinuity equation over the depth of the mixed- 
layer, assoeiated with the boundary eondition that the vertieal eomponent of the residual 
veloeity vanishes at the sea surfaee, leads to: 


dU ^ 
dx dy 


(7) 


For eonvenienee, in following equations, ^ will be written as . 


The next step eonsists of introdueing a eonvenient measure of the strength of the 
residual flow that enters (eseapes) the mixed-layer at a given sea surfaee density range. 
Consider a thin strip of fluid bounded by two nearby isopyenal surfaees (a and cr -i- Act ) 
as indieated by the shaded area in the sehematie in Fig. 4. The integral of over this 
area (5 ) is defined by: 

G{ct) = ^^w^JS ( 8 ) 

5 
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Figure 4. Schematic diagram illustrating the proposed diagnostic framework. 


The area of integration in Eq. (8) extends from one intersection with the coastal 
boundary to another if the buoyancy contours are blocked by the land [see the schematic 
diagram in Fig. 4], and it otherwise may include contributions from several ocean basins. 
For closed surface density contours, for instance in the re-entrant Antarctic Circumpolar 
Current, the integral pertains to the entire area between contours cr and cr -l- Act . 


The goal of the following derivation is to express G in terms of observed 


quantities. By using the relation between w^^^and U and V provided in Eq. (7), and 
Green’s theorem, the area integral G(cr)can be expressed by: 




dU ^ 
dx dy 


dS = j){Udy-Vdx) 


where dS is the boundary of the region S . Note that the direction for the loop integral is 
counterclockwise. 

Assuming no flux across the coastal boundaries for zonally blocked flows (or the 


exact cancellation of the zonal fluxes in case of a re-entrant flow), the integral in Eq. (9) 
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may be simplified to: 


G(cr) = J(t/dy -Vdx)- J {Udy-Vdx) 


( 10 ) 


(j+A<j 


where the integration along the isopycnals is in the positive v-direction. 

By introducing the unit vector normal to the buoyancy contours n = 

(10) can be conveniently rewritten: 

G(cr)= J U-ndl-^U-ndl 

(J+A(T a 


Vcr 

_ 

|Vcr„ 


Eq. 


( 11 ) 


with U =(U,V) and I being the arc length along the mixed-layer density contoursd 


Eq. (11) is readily interpreted as follows: the area-integrated vertical residual flux at the 
bottom of the mixed-layer in a steady state is equal to the sum of the lateral residual 
fluxes into this area. This statement is consistent with a simplified mass budget in the 
volume defined by the two isopycnal surfaces within the mixed-layer, as depicted by the 
schematic diagram in Eig. 5. This mass budget is here expressed as a simple volume 
budget since the density is almost constant within the considered volume and the 
contribution of evaporation/precipitation at the sea surface to the volume balance is 
neglected. 

Taking the limit Acr^0,Eq. (11) becomes: 

dF 

G(o-) = —-Ao- (12) 

^ ’ dcT 

where F{<j) = ^U-ndl (13) 


^ While the schematic in Fig. 4 and the related discussion represents the case in which density 
increases northward (the configuration relevant for the Northern hemisphere), it can be readily shown that 
our formulation applies to both Northern and Southern hemispheres. 


11 




Mixed-layer 

depth 


Figure 5. Physical interpretation of Eq. (11): mass budget in the volume delimited by the 
isopycnal surfaces cr and cr -i- Act within the mixed-layer. 


Eq. (6) may be written as t7 - n = -—-— and then substituted in Eq. (13) to obtain 


the following expression of F: 

F(a) = f 

^ \V a 

G \ m 


(14) 


Since the air-sea fluxes and the sea surface density are easily accessible from data, the 
preceding analysis offers a simple procedure for estimating the total residual volume flux 
that enters (escapes) the mixed layer from (into) the ocean interior within the surface area 
S(crj < cr < 0 - 2 ) that is bounded by the two density contours cr^ and cr^ (cTj > crj : 

y(cr;,cr2)= JJ w^^^dS = ^—^dcT=F{cT^)-F{cT,) (15) 


where F (cr) is given by Eq. (14). 
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B. ALTERNATIVE METHOD 

Instead of computing F(cr) along the ispopycnals cTj and cTj, y(cri, 0 - 2 )may be 

obtained by an alternative method using the area integration over the surface bounded by 
two isopycnals. This alternative technique is similar to that in Walin’s (1982) 
formulation. 


By introducing the expression for F [Eq. (14)] into Eq. (11), we evaluate the 
vertical residual flow that enters (escapes) the mixed-layer over the area S (defined above 
and sketched in Eig. 4): 


G(cr) = II= F(cr-I- Act)- F(cr) = | . 

S ij+A <7 \ ^m cr 


B, 


|V, 


'0 




(16) 


Alternatively, Eq. (16) may be expressed as: 


G(a) = 


Aa 


•’ Vcr„l J 






a. 


m\ J 


= ^f I B^Ldl-\B^Ul 


. (7+A(T 


(17) 


where L = ^-r is, approximately, the distance between the two isopycnal contours cr 

Vcr 

I ^ \ 

and cr -I- Act . Therefore, Eq. (17) becomes: 


G(a) 


1 

Act 


II Bg (cr-r Acr)d5 -|| Bg (cr)d5 

s s J 


(18) 


Taking the limit Act ^ 0, we arrive at: 




ea‘ 


Aa 


(19) 


and finally obtain the alternative expression forV (a^,a 2 ): 


7G(cr) 

\ ^ ’ da 

Act 

d 

( 

rr 


d 

r 


da 

JJ B,dS 

da 

JJ B,dS 



J 

(7=(T2 




( 20 ) 
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where 5 (cr) is the area bounded by the mixed-layer density contour <j„=<j and a 
reference isopycnal <Jq located in the same hemisphere and verifying cTq < cr. Note that 
the exact value of that reference density does not affect Eq. (21); therefore, for 
convenience, we may define cr^ as the minimum density with distinct non-intersecting 
northern and southern contours, as depicted in Fig. 6a. 

The alternative method leads to the following expression for F\ 



While this expression is identical to the water mass transformation function 
(Walin 1982), its role and interpretation in the model derived here are less conventional. 
First, we note that unlike Walin’s model and its previous extensions, this formulation 
takes into account the volume transfer driven by the adiabatic component of mesoscale 
eddies. More importantly, in this study, function F(cr) is defined [Eq. (13)] in terms of a 
horizontal volume flux in the mixed-layer crossing a given isopycnal surface <j. This 
connection makes it possible to calculate the net residual volume flux that escapes 
(enters) the mixed-layer from (into) the ocean interior within an area bounded by two 
density contours [Eq. (15) and Eq. (20)]. Walin (1982) did not separate the water mass 
transformation within the mixed-layer from that occurring in the ocean interior. 

It is also essential to note that the foregoing model makes it possible to distinguish 
the density forcing in the Northern and Southern Hemispheres. On the contrary, the 
original Walin’s formulation involved the total transformation at a given interior 
isopycnal surface, in which contributions from the air-sea density fluxes in both 
hemispheres were grouped together. 
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C. DIAPYCNAL VOLUME FLUX 

As defined above, cTq denotes the minimum value of mixed-layer density for 
which the studied basin contains two non-intersecting isopycnal contours. As shown in 
Fig. 6a, for every cr > cTq there are two distinct values of F corresponding to the southern 

and northern density contours; these are denoted as F^{(j) and Fj^(a) respectively. 



Figure 6. Schematic representation of the density structure in a zonally bounded ocean. 

(a) Mixed-layer density distribution. For each value of cr > cTq there are two non¬ 
intersecting contours of surface density; one in the northern part (indicated by the 
gray shading) and one in the southern part (not shaded), (b) The vertical 
meridional cross-section of density. The residual flux entering the density layer 
bounded by the isopycnals cTj and <J 2 from the mixed layer ( SV^ + ) is 

balanced by the diapycnal flux [ (a 2 ) - (cTj ) ] ■ 


Since the residual flow is assumed to be non-divergent, the total residual flux 
(14)that crosses the isopycnal cr in the ocean interior (see the schematic in Fig. 6b) is 
exactly equal to the flux from the ocean interior into the mixed-layer within the area 
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defined by the condition a^^<a : 

V^(a) = JJ G(a')da' = F^(a) + F^{<j) for cj>cj^. (22) 

<7'<a 

The sign convention in Fig. 6 is such that is positive (negative) for an upward 
(downward) integrated flux across the isopycnal in the ocean interior. 


D. ADIABATIC ADVECTION 

The next step (and the key strength of the proposed technique) is to quantify the 
magnitude and pattern of the isopycnal pole-to-pole MOC component, caused by the 
adiabatic interior advection, from the surface data. For this, the residual volume flux 
entering (escaping) the mixed-layer from (into) the interior is computed separately in the 
Northern and Southern Hemispheres. 

The highest value of potential density that can be found in the mixed-layer is 
denoted by and cr is an arbitrary value within the interval [crQ,crjjj^] . For each 

hemisphere, Eq. (15) is used to express, in terms of F , the net residual flux at the bottom 
of the mixed-layer within a density range of [cr, cr^^ ], as follows: 




jj WresdS = F^ (o-„^) - iv (cr) 

S(o-<o-’<o-„„) 

north 


5(o-<o-'<C7„„) 

south 


(23) 


Since the residual flow is non-divergent, the sum of all fluxes entering the interior 
volume bounded by the isopycnals cr^^ and cr is zero [see the schematic in Fig. 6b]: 


yA^^j-yA^)-yNicT)-v,{cj)=Q ( 24 ) 

Note that in the absence of any cross-isopycnal flows, Eq. (21) would reduce to 
V^(cr) =-y^(cr), which simply means that all the water subducted in the Southern 
(Northern) Hemisphere upwells in the Northern (Southern) one. In this case, the 
isopycnal flux throughout the density layer [cr,crj^^^] is uniform and the adiabatic volume 


transport is (cr) = (cr) = -V^ (cr). 
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However, in reality, the diapycnal flux across any isopycnal surface is not zero, 
and consequently V^(cr) is not equal to VjCcr). Therefore, the isopycnal flux cannot be 

precisely determined from and V^, and we need to define a way to assess the 
magnitude of V^. It is reasonable to assume that representative values of lie in 
between V^(cr) and -I4(cr) . Thus, without any additional assumptions, the simplest and 
most consistent estimate of the isopycnal flux 1/ within the density range [a, ] is 

given by the linear average of the two net residual fluxes at the bottom of the mixed- 
layer: 

( 25 ) 

Note that, with this definition, a positive value of is associated with a 

downward residual flux at the base of the mixed layer in the Southern Hemisphere and its 
corresponding upward flux in the Northern Hemisphere. While this definition seems 
somewhat arbitrary at first sight, its application to the global overturning circulation 
model, described in section IV-D, indicates that it is surprisingly accurate in predicting 
the actual inter-hemispheric volume exchange.^ 

Since the length of the mixed-layer density contours inevitably reduces to zero as 
the value of density approaches its highest possible value (o' ^ we may assume 

that ~ ~ 0 lu Eq. (23), and therefore Eq. (25) reduces to: 

r,(CT) = i[Fj(<T)-f„(<T)] (26) 

It is important to note that the expression above is most relevant for the case of a 
global ocean circulation. On the contrary, for computations within each basin, where we 
only consider the isopycnal contours zonally bounded by land, a better way to define 

and Vj^ than that expressed in Eq. (23) may be: 

2 Note that cr in the preceding formulation cannot be less than , and therefore Eq. (25) does not 
describe the overturning in the upper tropical thermocline where a < a^. However the region 
corresponding to this density range occupies only a small fraction of the thermocline, with a maximum 
depth of 200 m, and thus its relative contribution to the overall overturning is expected to be negligible. 
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(^) = jj ^resdS = iV (^) - (^^0 ) 

5{a-o<cr'<cr) 
north 

(a) = jj =F, (a) - (a^ 

5((Tq<(j'<(j) 
south 


With this definition, the integration is carried out over the area defined by the 



mixed-layer density range [cro,cr], which is more representative of the considered 


oceanic basin than the density interval [cr,cr^^,.]. Moreover, with the alternative 


definition of and provided in Eq. (27), the adiabatic flux defined in Eq. (25) 
now represents the isopycnal volume transport within the density range [crg,cr]. 


E. ADIABATIC FLUX RATIO 

(cr) in Eq. (26) represents the integrated isopycnal volume flux over the 
density interval [cr, cr^^^ ]. Therefore, the fraction of this adiabatic volume flux within the 

dV 

density layer [cr,cr-i- Acrl is —- Acr, where Act is the density increment. 

5cr 

One of the objectives of this study is to compare the diapycnal and isopycnal 
components of the overturning circulation, and some quantification of their relative 
strengths is obtained through the following parameter: 


/ = 




(28) 



= rms 


\d(7/ 

^0«^<^m3x 

d(j 


rms 

(7(i«7«7^. 


Iy,(a) I, and Aa=a-(7.,. 


Parameter / is referred to hereafter as the adiabatic flux ratio. If /»l, the 
interior flow is largely adiabatic, whereas / « 1 indicates the importance of the interior 
diffusive processes for maintenance of the MOC. 
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F. EKMAN FLUX 

As defined above, the residual flux includes the mean circulation and the eddy- 
contribution. The Ekman flux, which is generated by the mechanical forcing by the 
surface winds, dominates the lateral mean flux in the mixed-layer. Therefore, in regions 
where the eddy activity is weak, the residual flux may be almost equal to the Ekman flux. 


In particular, the Ekman pumping velocity at the bottom of the Ekman layer is 
given by: 


A 

dx 


( T ^ 

‘'Qy 

d 

fr ^ 

‘-Ox 

= curl 

f - 3 

T 

\PJ ) 

dy 

yPof J 


V Pof y 


(29) 


where and Tq^, are respectively the zonal and meridional components of the wind stress 
f at the sea surface, p^^ is the reference seawater density, and/is the Coriolis parameter. 


In a similar way to the expression of the net residual flux at the bottom of the 
mixed-layer [Eq. (23)], we may define for each hemisphere the net Ekman flux at the 
bottom of the Ekman layer within a density range of [cr, ], as follows: 


jj w^dS 

S(o-<o-'<o-„„) 

north 


< 

~ JJ w^dS 

S(cr<o-'«T„„) 

south 


(30) 


As mentioned in Section II-D, in the case of a specific oceanic basin, where only 
the isopycnal contours zonally bounded by the land are considered, we estimate the 
overturning circulation in the density range [crQ,cr]. Therefore, the alternative 
expressions of the northern and southern Ekman fluxes are in this case: 

^N-EK (‘^) ^ JJ W^dS 




north 


JJ 


( 31 ) 


south 
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Finally, by analogy to the previous definitions relative to the residual flux, we 
may define the Ekman adiabatic advection (wind-driven component of 1/ ): 

= (32) 

The comparison of and V^(a) allows us to assess the role of the wind 

stress in the establishment and maintenance of the MOC. 
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III. PARAMETER DEEINITION AND DATABASE DESCRIPTION 


The technique described in the previous chapter offers a quantitative description 
of the two components (adiabatic and diapycnal) of the MOC based on sea surface 
observations. The next step consists of choosing the datasets in order to provide an 
accurate assessment of the oceanic overturning. 

Since the MOC varies on time-scales which greatly exceed the typical scale of the 
air-sea fluxes (diurnal cycles), the model used in this study must rely on long-term means 
of observed parameters, and the computations carried out below use climatology datasets. 
This approach is consistent with the variability of the studied phenomena and secondly it 
allows the analysis of climate variations, which is one of the main motivations for 
studying the MOC. 

The computations in this study use annual means for all quantities and it is 
essential, throughout the data description, to examine whether the climatology averages 
accurately represent the reality. 

A. MIXED-LAYER DENSITY 

I. Definition 

One of the quantities used in the computation of the residual flux is the mixed- 
layer density anomaly (cr^). The mixed-layer density is conceptually well-defined: the 

intense mixing within the mixed-layer (ML) makes the seawater density vertically 
uniform. Therefore, in the conceptual models, the density field at any depth above the 
pycnocline accurately represents the ML density distribution. 

However, in reality, the annual-mean density is not vertically uniform in the 
mixed-layer and the choice of one particular depth to characterize the ML properties is 
not trivial. This choice is made more difficult by the large seasonal and spatial variability 
of the ML depth. Indeed, our theory uses a spatially uniform model of the mixed-layer 
depth for the computation of the density gradient that is needed in the calculation of the 
MOC components [Eq. (14)]. 
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Fig. 7 shows the spatial distribution of the ML depth: the contours represent the 
annual means provided by the Levitus 94 climatology. This map indicates the large 
difference between the subtropical and polar regions. These differences are significantly 
reduced here because the plotted data are the annual averages and this filters out the 
seasonal variations. The observations in Fig. 7 suggest that the mixed-layer density can 
be approximated by the density at the sea surface. This idea is supported by the plot in 
Fig. 8, which gives the annually- and zonally-averaged ML depth with its temporal 
standard deviation computed over the monthly mean values. 

On the other hand, as shown in Table 1, the temporal variations of density tend to 
significantly increase as the chosen depth is closer to the sea surface. Thus, in order for 
cr(/z) to accurately represent the mixed-layer density distribution throughout the year, 

we have to choose depth h large enough to have a low temporal standard deviation. 
Moreover, a deeper level within the ML is more relevant for describing the density layer 
of the MOC in the ocean interior. 

Thus, a compromise must be found in order to select a depth which is almost 
always and everywhere within the mixed-layer and at which the density temporal 
variability is relatively low. The depth chosen for this study is 40 meters, based on Fig. 8 
and Table 1. Note also that the 40 m depth is not within the ML only in the equatorial 
region and in the Arctic Ocean. The equatorial area is excluded from our analysis: we 
only consider regions where o (Jq and the technique derived in the previous chapter 
does not describe the overturning in the upper tropical thermocline where cr < cTq . Since 

the North Pole region is mostly covered by ice, the density forcing at the sea surface there 
will not be taken into account in our computations. 
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Figure 7. Annual mean of the oceanic mixed-layer depth given in meters and based on 
Levitus 94 climatology. White areas represent the continents and coastal regions 
where the mixed-layer depth is less than 20 m. 


In order to more accurately represent the ML density, the small effect of the 
pressure on the density must be removed; hence, the potential density anomaly at 40 m is 
used for the computations. The potential density anomaly (cr^) is related to the potential 

density (p^in kg trT ^) as follows: 

Pg =1000 -I- CTg 

Note that using the potential density field to represent the ML density brings our analysis 
closer to the conceptual models of the mixed-layer, since this parameter has a lower 
vertical gradient than density. 
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Figure 8. Annual mean and temporal standard deviation of the zonally-averaged mixed- 
layer depth. The plotted data are based on Levitus 94 climatology. The solid line 
indicates the zonally- and annually-averaged ML depth distribution and the dotted 
lines point out its temporal standard deviation computed from the monthly ML 
depth means. The dashed horizontal line highlights the 40 m depth. 


Depth (m) 

Standard deviation of potential density 

(kgm''') 

0 (surface) 

4.7945 

10 

0.3995 

20 

0.3565 

30 

0.3184 

40 

0.2670 

50 

0.2403 

75 

0.1720 

100 

0.1349 


Table 1. Mean temporal standard deviation of potential density at each characteristic depth 
of the mixed-layer. These statistical values are computed from the Levitus 94 
monthly means of potential density at all grid points of the whole ocean. 
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2. Levitus Climatology 

The oceanic database used in this study, for the potential density distribution at 40 
m, as well as temperature and salinity needed in the calculations of the density flux at the 
sea surface, is Levitus 94 climatology. This database consists of a world atlas of 
objectively analyzed mean fields of major oceanic parameters at the annual, seasonal, and 
monthly time scales. 

All the analyzed fields are offered on a one-degree latitude-longitude grid at 
selected standard depth levels from the sea surface to 5500 m depth. The standard levels 
relative to the ML are 0, 10, 20, 30, 50, 75, 100, 125, 150, 200, 250, and 300 meters. 
Therefore, for the 40 m potential density anomaly distribution, we have linearly 
interpolated the 30 m and 50 m analyzed fields. 

The atlas datasets have been produced by the Ocean Climate Laboratory (OCL) at 
the National Oceanographic Data Center (NODC), which is supported by the National 
Oceanic and Atmospheric Administration (NOAA) Climate and Global Change program 
to provide quality controlled oceanographic databases. In particular, they include an 
objective analysis of all the in situ measurements for the period 1900-1992. 

The Levitus climatology is considered by all oceanographers to be one of the 
most reliable databases for the oceanic parameters. It is accessible through the website 
http://iridl.ldeo.columbia.edU/SOURCES/.LEVITUS94 in NetCDE (Network Common 
Data Eorm) format. 

3. Potential Density Field 

The annual mean distribution of potential density anomaly at 40 meters 
{(Tg (40m)) from Eevitus 94 climatology is given in Eig. 9. The isopycnal contours <Jq in 
both hemispheres are highlighted on the map. Based on the definition of given in the 
previous chapter and the monthly mean fields of cr^ (40m) [not shown here], the value of 
cTq has been estimated at cr^ = 24.2. This value corresponds to the minimum potential 
density with distinct non-intersecting northern and southern contours over the whole 
ocean and throughout the year (i.e., on every chart of cr^ (40m) monthly means). 
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Figure 9. Annual mean distribution of potential density anomaly at 40 m depth. The 
dimensionless displayed data have been computed from Levitus 94 climatology. 
The thick solid lines point out the northern and southern contours of 

(cTq = 24.2). Dotted lines indicate partitioning of the world ocean into the 
individual basins. 


B. 


DENSITY FLUX AT THE SEA SURFACE 


1. Definition 

The second surface quantity used in the computation of the residual flux is the air- 

sea density flux (B in kg ). As shown previously in Eq. (1), this flux is composed 

of two terms, the thermal part (5j) and the freshwater component (fij): 

^ aH _ {E-P)S 
B = - + ^ 




aH 


1-5 


C„ 


^2 = M 


{E-P)S 

1-5 


(33) 

(34) 
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Note that many studies dealing with fluxes at the sea surface refer to the buoyancy 
flux: this latter may be obtained from the density flux by simply multiplying B by the 

constant —— , with g being the gravity constant (9.8 m/ 5 ^) and the standard 
Pq 

seawater density (l025 kg . Indeed, buoyancy b is defined as b = — —p , where p 

Po 

is the seawater density. 

The parameters , a and , which are, respectively, the specific heat capacity, 

the thermal and saline expansion/contraction coefficients of seawater, are not constant 
and depend on temperature, salinity, and pressure: 

a = a[T,S,p) 

■fi = p(T.S.p) 

c„ = c,(r, 5 ,p) 

Therefore, the temperature (7) and salinity (5) distributions at the sea surface 
(where the hydrostatic pressure p is equal to 0), which are available in Levitus 94 
climatology, will be used for determining the air-ocean density flux. The three 
coefficients are then calculated on the basis of the UNESCO algorithms for computation 
of fundamental properties of seawater (1983). 

The heat flux from the atmosphere into the ocean (H), the evaporation (E) and 
precipitation (P) rates are obtainable in meteorological databases. In particular, the 
following fields may be found in the different climatological databases: 

Net solar energy flux at the sea surface (2^); 

Net longwave radiation at the sea surface (2 l ); 

Latent heat flux (7/^); 

Sensible heat flux (7/^); 

Total precipitation at the surface (7*), which is sometimes divided into two 
fields (large-scale and convective precipitation). 

Note that, for the four thermal fluxes , the convention used here 

associates positive values to upward fluxes (from the ocean into the atmosphere). 
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The five parameters described above allow us to compute the heat flux from the 
atmosphere into the ocean {H in W wT^) and the evaporation and precipitation difference 
(E-P in m ^ ‘) field through the following equations: 

(35) 

and E-P = ^ — P (36) 

where is the specific latent heat for evaporation w 2.5x10^ J ^. 

2. Atmospheric Databases 

Four different databases have been used to compute H and E-P. We compare their 
climatology and the corresponding predictions of the MOC components. This comparison 
answers two essential questions: the quality control of the results and their sensitivity to 
the observational fields. 

• The European Center for Medium-range Weather Forecast (ECMWF) Re- 
Analysis (ERA-40) is a global analysis of many observations and satellite data 
streams for the period September 1957 through August 2002. This dataset 
contains 56 surface and single-level variables on an equally-spaced global 2.5 
degree latitude-longitude grid, with all variables reported four times a day. In 
particular, the ERA-40 atlas, the web-version of which was released on 6 
February 2006, describes the climate from 1979 through 2001, which is the period 
with the best and most time-consistent product quality for the entire globe. For 
this reason, this study uses monthly means of each parameter listed above from 
January 1979 through December 2001. The data may be downloaded from 
http://dss.ucar.edU/datasets/dsll8.0 . Note that the data are provided there in GRIB 
format (binary code for exchange of processed data). 

• The joint National Centers for Environmental Prediction (NCEP) and 
National Center for Atmospheric Research (NCAR) re-analysis project has 
performed data assimilation using historical data from January 1948 to the 
present. Particularly, the NCEP-NCAR Climate Data Assimilation System 1 
(NCEP-NCAR CDAS-1) provides monthly means of each needed parameter from 
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January 1949 up to the current month on an equally-spaced global latitude- 
longitude grid containing 94x192 points. In order to be consistent with the 
ECMWF climatology, we use the subsets from January 1979 through December 
2001 in this study. These datasets may be obtained at 
http://iridl.ldeo.columbia.edU/SOURCES/.NOAA/.NCEP-NCAR/.CDAS-l in 

NetCDE format (Network Common Data Form). 

• The Global Air-Sea heat and momentum flux Climatology from the 
Southampton Oceanography Centre (SOC-GASC97) provides climatologic 
sea surface fluxes based on in situ reports within the CO ADS lA (Comprehensive 
Ocean Atmosphere Dataset lA), a global dataset containing of order 30 mi llion 
surface observations from ships an buoys collected over the period 1980-1993. It 
is important to note that the quality of the fields is known for having a strong 
spatial dependence that reflects the global distribution of ship observations. 
Quality is likely to be high in the well-sampled North Atlantic and North Pacific 
oceans but, it decreases in the Southern Hemisphere. In particular, south of 40 S, 
the errors in the fields are likely to be large and the existence of spurious features, 
which have been generated during the objective analysis of the original raw 
distributions, must be pointed out. The data are available on an equally-spaced 
global one-degree latitude-longitude grid, in NetCDE format, at 
http://iridl.ldeo.columbia.edU/SOURCES/.SOC/.GASC97 . The dataset 
documentation is accessible from this website. 

• The Da Silva atlas of Surface Marine Data 1994 (SMD94) provides the 
monthly climatology of the near-surface atmospheric properties over the global 
ocean. Covering the period of January 1945 through December 1989, this 
database is made from individual observations in COADS lA and an objective 
analysis using Barnes’ weight function, which is the same response function used 
in Eevitus 94 Climatology. It was produced at the University of Wisconsin- 
Milwaukee (UWM) by A. M. da Silva and C. C. Young in collaboration with S. 
Eevitus from the National Oceanic and Atmospheric Administration (NOAA) / 
National Oceanographic Data Center (NODC). The datasets are available on an 
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equally-spaced global one-degree latitude-longitude grid, in NetCDF format, at 
http://iridl.ldeo.columbia.edU/SOURCES/.DASILVA/.SMD94 . The dataset 
documentation is accessible from this website. 


3. Distribution of the Density Flux at the Sea Surface 

Based on Eqs. (1), (33), (34), (35), and (36), the density flux at the sea surface has 
been computed from each of the four databases. Eigs. 10, 11, 12, and 13 show the density 
flux distributions respectively from ECMWE re-analysis, NCEP/NCAR re-analysis, 
SOC-GASC97, and DaSilva-SMD94 climatology. In each figure, the thermal and 
freshwater components (respectively and ) are also displayed. 


Over most of the oceans, the density flux seems to be driven by thermal forcing. 
Eigs. 10, 11, 12, and 13 show that the distributions of total density flux B and its thermal 
part look similar. Moreover, by looking at their magnitudes, we see that is 

generally the dominant part of 5. This observation is supported by the zonal means of 
annually-averaged values that are given, respectively for each meteorological database, in 
Eigs. 14, 15, 16, and 17. 


Thus, from a global prospective, thermal forcing dominates the freshwater forcing 
and controls the density flux at the sea surface. However, as shown in Eigs. 18, 19, 20, 


and 21, which provide the annually-averaged distribution of log 


/ 

5 , 

\ 

V 

B2 



for each database. 


the zonal variability is large. In some places, especially in subtropical and polar regions, 
^2 greatly exceeds B^. These observations become more critical when we consider 
seasonal fields [not shown here]. 


Eigs. 10, 11, 12, and 13 show the same regions of high density flux: these are in 
general the Western Boundary Currents (WBC) and the equatorial band. Indeed, in the 
regions of the warm Gulf Stream, Kuroshio Current, Agulhas Current, and East 
Australian Current, the heat loss by the ocean into the atmosphere is large, which causes 
the upper ocean to gain density. These currents advect heat from the tropics towards the 
poles, and energy exchanges with the atmosphere occur in their paths. 
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Figure 10. Sea surface density flux from ECMWF climatology. Positive values correspond to 
density gain by the ocean. The displayed values are in kg (a) Total density 

flux B (b) Thermal component (c) Freshwater component B^. 
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Figure 11. Same as Fig. 10, but from NCEP/NCAR re-analysis (same period as ECMWF 
climatology: 1979-2001). 
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Figure 12. Same as Fig. 10, but from SOC-GASC97 climatology (period: 1980-1993). 
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Figure 13. Same as Fig. 10, but from DaSilva-SMD94 climatology (period: 1945-1989). 
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Along the equator, the upwelling induced by the Ekman divergence brings cold 
water to the upper ocean and leads to vigorous downward heat fluxes from the warm 
tropical air above the surface. The negative oceanic density flux, caused by thermal 
forcing, is enhanced by the large precipitation within the Inter-Tropical Convergence 
Zone (ITCZ) band, which is, on average, centered on the equator. 

Concerning the freshwater contribution, the large evaporation and the weak 
precipitation in subtropical regions generates a strongly positive contribution to the 
density flux at the surface. These regions are characterized by high sea-level pressure and 
persistent coverage of low-level non-precipitating clouds. On the other hand, over the 
equatorial band (mean position of the ITCZ) and in the sub-polar areas, precipitation is 
greater than evaporation, resulting in negative oceanic density flux. 

The four databases considered in this study provide the same overall picture of the 
surface density flux distribution. However, significant quantitative differences are noted 
in some regions for the density flux components and . 

It is generally agreed that ECMWF and NCEP/NCAR are the two most reliable 
among these four climatologic databases. A number of studies have compared these two 
re-analysis datasets, and many of them (for instance, Bromwich and Eogt 2004) have 
suggested that ERA-40 seems to have better skills than NCEP/NCAR regarding analyzed 
fields after 1978 (the beginning of the modern satellite era). This may be explained by a 
better assimilation of satellite data within ECMWF than in NCEP/NCAR, which is 
largely constrained by the station observational network. Thus, ECMWF is particularly 
useful for the high latitudes. 
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Figure 14. Zonal means of annually-averaged density flux B at the sea surface and of its two 
components andfij from ECMWF climatology (1979-2001). 



Latitude C^) 

Figure 15. Same as Fig. 14, but from NCEP/NCAR re-analysis (1979-2001). 
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Figure 16. Same as Fig. 14, but from SOC-GASC97 climatology (period: 1980-1993). 
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Based on the foregoing considerations, ECMWF re-analysis has been chosen as 
the reference meteorological database in this study. This choice is supported by the charts 
in Fig. 22, which display the differences between the annually-averaged density fluxes 
obtained from ECMWF re-analysis and those computed from the other datasets. Indeed, 
ECMWF appears to be a good compromise between NCEP/NCAR, DaSilva-SMD94, and 
SOC-GASC97. NCEP/NCAR provides higher density flux than ECMWF by the 
subtropical west coast of major continents and over a large part of the subtropical North 
Atlantic Ocean. On the other hand, DaSilva-SMD94 and SOC-GASC97 give smaller 
values than ECMWF over all tropical and subtropical oceans expect for some west coast 
regions. These features are observed for almost all latitudes, as indicated by the zonal 
means of these density flux differences in Fig. 23. 

Table 2 provides the global mean value of the density flux differences to assess 
the global accuracy of each meteorological database versus ECMWF. These values 
support the choice of ECMWF as the reference climatologic database for this study. 



Mean density flux difference 

(kgm-"s-‘) 

B (NCEP / NCAR )-B{ ECMWF ) 

1.8716x10 " 

B[SOC-GASC9l)-B[ ECMWF ) 

-1.3075x10 ® 

B (DaSilva )-B{ ECMWF) 

-1.3756x10 ® 


Table 2. Global averages of density flux differences between ECMWF re-analysis and the 
other databases. 
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Figure 19. Same as Fig. 18, but from NCEP/NCAR re-analysis (1979-2001). 
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Figure 20. Same as Fig. 18, but from SOC-GASC97 climatology (period: 1980-1993). 



Figure 21. Same as Fig. 18, but from DaSilva-SMD94 climatology (period: 1945-1989). 
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Figure 22. Differences in the surface density flux (a) ^B[NCEP/NCAR)-B^ECMWpyj 
(b) [b{S0C-GASC91)-B{ECMWE)] (c) [b{D aSilva)-B{ECMWE)]. 
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Latitude 


Figure 23. Zonal means of the differences in the surface density flux: 

[^B[NCEPINCAR)-B[ECMWF)\, [^B[SOC-GASC91)-B[ECMWE)'\, and 

[^B[DaSilva)-B[ECMWE)'\. 

C. MOMENTUM FLUX AT THE SEA SURFACE 

As described in section II-F, the computation of the Ekman flux requires the zonal 
and meridional components of the wind stress f at the sea surface, denoted by and 

respectively. Three databases have been used for the wind stress components: 

ECMWF re-analysis, NCEP/NCAR re-analysis, and DaSilva-SMD94 climatology, with 
the same periods of record as those defined for the density flux calculations. The 
computed Ekman flux provides a convenient check of the residual flux in regions where 
the eddy activity is weak. 
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IV. RESULTS AND ANALYSIS 


A. GLOBAL OCEAN 

Fig. 24 plots the zonally-averaged meridional cross-section, in which the 
climatological potential density anomaly contours, computed from Levitus94 database, 
are displayed. This cross-section is similar to the meridional cross-sections of individual 
oceans and shows a density distribution that is analogous to the schematic in Fig. 6b. 



Figure 24. Zonally-averaged potential density anomaly (dimensionless) for the upper 1000 
meters of the ocean. The data have been computed from the Levitus 94 
climatology. 


An advantage of applying the diagnostic model to the global configuration is that, 
unlike the individual basins, the water mass budget of the global ocean does not contain 
the contribution from the external sources or sinks. The density range considered for the 
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full ocean-wide residual flux budget extends from cr^ = 24.2, determined in the previous 


chapter, to = 28.3, the highest density of the ice-free sea surface area. 


1. Diapycnal Volume Flux 

The diapycnal interior flux (V^), computed from the four meteorological 
databases using Eq. (21) and Eq. (22), is plotted in Eig. 25. Its general pattern and 
average value of 3.5 Sv (one Sverdrup is defined by 1 Sv = 10® '), obtained from the 

ECMWE climatology (reference database), are consistent with the earlier estimates by the 
Walin-type calculations (Donners et al. 2005). 

Eor cr < 26.1, (cr) is positive, which implies an upward integrated flux across 

the (j isopycnal surfaces in the ocean interior. This result is a consequence of the strong 
low-latitude upwelling and the net upward air-sea heat flux in mid-latitude regions. 
Indeed, as indicated in Eq. (14), a positive density flux into the ocean (induced by a net 
upward air-sea flux) leads to a positive value of E(cr) in each hemisphere, which implies 
water mass transformation to higher densities within the mixed-layer. 

The diapycnal flux then changes sign and remains negative over the density 
anomaly interval [26.1, 27.6]. This may be interpreted as a convergence of the residual 
flow within the mixed-layer in this density range. The negative values of are 

associated with the strong sinking motion of Antarctic Intermediate Waters (AAIW) and 
Sub-Antarctic Mode Waters (SAMW), which occurs along the Antarctic Circumpolar 
Current (ACC), characterized by a potential density anomaly of about 27; indeed, the 
cross-isopycnal flux has a negative peak value at cr« 27.1. 

Eor cr > 27.1, increases with cr and becomes positive from cr = 27.6 . In fact, 

the strong downwelling along the ACC is compensated by the large-scale upwelling of 
North Atlantic Deep Water (NADW) and North Pacific Deep Water (NPDW), which 
takes places in the Southern Ocean, between the ACC and the Antarctic continent. Thus, 
our diagnostic results in Pig.25 are consistent with the distribution of water masses over 
the whole ocean, sketched in Pig. 2. 
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Figure 25. Diapycnal interior flux in the global ocean computed from the four 
meteorological databases. 


Fig. 25 also highlights the strong discrepancies among the I 4 values obtained 
from the four climatologic databases. Only ECMWF and NCEP/NCAR datasets result in 
profiles consistent with the previous results (Donners et al. 2005). However, 

NCEP/NCAR values are comparatively too high at low density anomalies (cr < 26). 

SOC-GASC97 and Da Silva SMD94 databases generally provide inconsistent 
y^ profiles. This result could be expected, since, as described in the previous chapter, the 
air-sea fluxes in these two databases are questionable for a large part of the whole ocean. 

Table 3 gives the average cross-isopycnal flow differences between values 
obtained from ECMWE re-analysis and those computed from the other databases, over 
the potential density anomaly interval [24.2, 28.3]. This table should be compared with 
Table 2, which provides the global averages of density flux differences between ECMWE 
climatology and the other databases. The comparison shows that the diapycnal volume 
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flux estimates are critically dependent on the density flux accuracy. Indeed, as shown in 
Fig. 26, the absolute diapycnal volume flux accuracy is empirically related to the absolute 
mean density flux accuracy, by the following equation: 

V^dijf = 5.38x10"' Bdijf ' - 7.48x10“" Bdiff ' +3.90x10"' \Bdijf\ (37) 

where V^diff is the average diapycnal volume flux accuracy, Bdijf is the mean density 

flux accuracy, and all the coefficients in Eq. (37) are measured in the corresponding MKS 
units. 



Mean diapycnal volume flux 
difi'erence (m's 'j 

(NCEP / NCAR ) - ( ECMWF) 

5.0470x10' (~5 5v) 

V^{SOC- GASC91 ) - y^ ( ECMWF ) 

-4.3530x10" (~-43.5 5v) 

y^ ( DaSilva ) - y^ ( ECMWF ) 

-5.2320x10" {--51 Sv) 


Table 3. Means of diapycnal volume flux differences between ECMWF re-analysis and the 
other databases. 
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Figure 26. Empirical relation between the absolute mean density flux accuracy and the 
absolute average diapycnal volume flux accuracy. The red dots represent the 
differences in the density flux means and in the associated diapycnal volume flux 
averages between, from left to right, NCEP/NCAR and ECMWE, SOC-GASC97 
and ECMWE, and Da Silva SMD94 and ECMWE, respectively. 


2. Adiabatic Advection 

Eig. 27 gives the distributions of the along-isopycnal interior flux (E^), computed 
from the four meteorological databases using Eqs. (21), (23), and (25). The plots display 
the residual (total) fluxes and their Ekman components. These plots should be interpreted 
with caution, since they represent the integrated isopycnal interior flux over the potential 
density anomaly range . Thus, changes in the direction of the MOC circulation 

(northward or southward) are observable by maxima and minima in the profiles. 


The direction of the overturning at each cr is determined by sign of 


d(j 


Since is defined by V^(a) = ^[V^(cr) -^^(cr)] [Eq. (25)], an increasing (decreasing) 


from down to ct^ corresponds to a northward (southward) flux. 
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Therefore, the total pattern makes it simple to identify the individual 
components of the isopycnal MOC. Most notably, the southward flux of the North 
Atlantic Deep Water (NADW) is reflected in the strong decrease in from about 

<7 = 21.6 to cr = 27.1. The northward moving water masses originating in the ACC 
(Antarctic Intermediate Water and the Sub-Antarctic Mode Water) also have a clear 
signature through the significant increase in 17 from about cr = 27.1 to cr = 26.1. 

Based on the residual (total) isopycnal flux profile computed from ECMWF 
climatology, the sketch in Fig. 28 summarizes the mean adiabatic MOC circulation for 
the whole ocean. Fig. 29 is a schematic of the inferred zonally-averaged meridional 
cross-section of the whole ocean. 
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Figure 27. Adiabatic interior flux \7 in the global ocean computed from the four 
meteorological databases. 
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The profiles in Fig. 27 shows the pattern of the meridional overturning that is 

generally consistent, in terms of magnitude and direction, with earlier numerical and 
observational studies. 

It is also important to point out that signs and magnitudes of 1/ computed from 
the total flux and their wind-driven components based on the Ekman transport are poorly 
correlated. One of the interpretations of this result is that the global MOC adiabatic 
component may not be mainly driven by mechanical forcing generated by winds, and the 
eddy contribution is essential for the maintenance of meridional overturning. 



Figure 28. Schematic of the adiabatic meridional overturning circulation in the global ocean. 


The 1/ profile computed form the reference climatology (ECMWF) indicates that 

the integrated adiabatic advection over the whole density range is almost zero. 

This implies that, as expected, the northward and southward interior fluxes balance each 
other in this density interval, lending more credence to the proposed diagnostics. 
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Figure 29. Schematic of the zonally-averaged meridional cross-section of the global ocean 
and the associated MOC. Red arrows label southward adiabatic advection and 
blue arrows depict northward adiabatic advection. 


In addition, Fig. 27 highlights the discrepancies among the values obtained 
from the four climatologic databases. The differences in the sea surface density flux 
between these datasets do not result in as large changes in 1/ profile as those in the 
computations. Moreover, we still observe that, similar to I 4 plots (Fig.25), the profiles 

based on ECMWF and NCEP/NCAR climatology are relatively close to each other, while 
those calculated from SOC-GASC97 and Da Silva SMD94 have significant differences. 

Table 4 gives the average along-isopycnal flux differences between values 
obtained from ECMWE re-analysis and those computed from the other databases, over 
the potential density anomaly interval [24.2, 28.3]. As in section IV-A-1, results in this 
table should be compared with those in Table 2. This comparison shows that the 
estimates of isopycnal volume flux are much less dependent on the density flux accuracy 
than are those of diapycnal volume flux. Indeed, the differences among the different 

databases are much lower that the corresponding differences. Unlike , no coherent 
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empirical function relating the absolute mean density flux accuracy and the absolute 
average adiabatic flux accuracy was found. This result could be expected since signs of 
mean density flux differences and mean isopycnal volume flux differences are not 
correlated [see Tables 2 and 4]. 



Mean isopycnal volume flux difference 
(mV) 

(NCEP / NCAR ) - ( ECMWF) 

-1.2506x10' (~-lSv) 

V^{SOC- GASC91 ) - y^ ( ECMWF ) 

-4.0233x10' (~-4Sv) 

y^ (DaSilva )-V^{ ECMWF ) 

-1.1483x10" (~-nSv) 


Table 4. Means of isopycnal volume flux differences between ECMWF re-analysis and the 
other databases. 


3. Adiabatic Flux Ratio 

Table 5 gives the adiabatic flux ratios computed from each meteorological 
climatology using Eq. (28). Here, the adiabatic flux ratios from ECMWE and 
NCEP/NCAR are about 3.5 and 3.0, respectively. Their magnitudes support the theory 
that the global Meriodional Overturning Circulation has a dominant adiabatic component. 
The SOC-GASC97 and Da Silva SMD94 values are provided here only for references 
without further interpretation, since their corresponding and 1/ profiles appear to be 
inconsistent with the real oceanic circulation. 


Database 

Adiabatic flux ratio / 
(dimensionless) 

ECMWF 

3.5018 

NCEP/NCAR 

3.0288 

SOC-GASC97 

1.4552 

Da Silva SMD94 

2.3247 


Table 5. Adiabatic flux ratios for the global ocean computed from the four meteorological 
databases. 
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B. ATLANTIC OCEAN 

If we limit our study of the Atlantic Ocean to the density range for which all 
isopycnals are zonally bounded by land throughout the year, the potential density 
anomaly interval reduces to [24.2, 25.4]. These values are obtained from the analysis of 
the monthly mean distribution of potential density [not shown]. As mentioned previously, 
confining the analysis to this density range permits an accurate description of the upper- 
ocean meridional circulation. This assures the absence of major external sources of water 
as assumed in the model formulation (section II). 

V'^(cr) and over the sigma range [24.2, 28.3] (used for the global ocean 

computations) have also been computed in order to assess the Atlantic contribution to the 
global MOC, although their pattern in the range [26.5, 28] could be affected by the 
contribution from the Antarctic Circumpolar Current (ACC). 

1. Diapycnal Volume Flux 

The diapycnal interior flux (V^) in the upper Atlantic Ocean, computed from the 

four meteorological databases using Eq. (21) and Eq. (22), is plotted in Eig. 30. Its 
average over the potential density anomaly range [24.2, 25.4] is approximately 6.7 Sv 
from the ECMWE climatology (reference database). 

Eig. 30 shows that is mostly weak over the density interval [24.20, 25.16]: its 

values remain between -9 Sv and -i-8 Sv (ECMWE). The diabatic volume flux abruptly 
increases to reach its maximum magnitude (37 Sv) at aboutcr = 25.34. This peak value 
corresponds to the high upward net heat flux at the surface of the Gulf Stream. Indeed, 
the mean potential density anomaly at 40 m depth within the Gulf Stream is cr » 25.35 . 


52 




Figure 30. Diapycnal interior flux in the Atlantic Ocean computed from the four 
meteorological databases. 


Similar to the global diapycnal volume flux distribution, we note large 
discrepancies between the estimates from different databases. Results from ECMWF 
and NCEP/NCAR re-analyses show many similarities, but NCEP/NCAR density forcing 
seems to be too high with regard to values. Concerning SOC-GASC97 and Da Silva 

SMD94 climatologies, the results are again unrealistic, which we attribute to the lack of 
observational data for the Southern Atlantic Ocean. This last point is illustrated in Pig. 
31, where the values of F{(j) = ■ ii dl [defined in Eq. (13)] in the Northern and 

<J 

Southern Hemispheres are given for the four climatologic databases: the differences in 
F^ values are much lower than those in F^ estimates. 
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2. Adiabatic Advection 

Fig. 32 gives the along-isopycnal interior fluxes {VJ, computed from Eqs. (21), 
(25), and (27) using the four meteorological databases. Here, (cr) estimates represent 
the integrated isopycnal interior flux over the potential density anomaly range[crQ,cr] 
instead of over the density interval [cr, cr^^ ] (as in the global ocean computations). 


Changes in the direction of the MOC circulation (northward or southward) are 
still observable by maxima and minima in the profiles, but, in this case, the direction 

dV dV 

of the adiabatic advection at each cr is given by the sign of —-: —- > 0 means a 

5cr d<j 


northward flux, and < 0 corresponds to a southward flux. (cr), calculated from 

ECMWF (reference climatology), shows a weak southward integrated isopycnal flow for 
24.20 < cr < 25.17 and then a strong northward flow for 25.17 < cr < 25.38 . This points 


out the predominantly northward residual flow in the Atlantic upper thermocline. Figure 
32 also highlights the strong effect of the large and positive density flux (into the ocean) 
in the Gulf Stream region on the adiabatic flow in the upper Atlantic Ocean. 


Concerning the Ekman fluxes, the plots suggest that the MOC in the Atlantic 
Ocean is not driven by the local mechanical forcing (wind stress). This result is similar to 
that obtained for the global thermohaline circulation. 
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Figure 31. 


Cross-isopycnal volume flux within the mixed-layer, (a) In the Northern 
Hemisphere (cr)^. (b) In the Southern Hemisphere (*^)) • 
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Figure 32. Adiabatic interior flux in the upper Atlantic Ocean computed from the four 
meteorological databases. 


3. Adiabatic Flux Ratio 

Table 6 gives the adiabatic flux ratios for the Atlantic basin. This parameter is 
relatively high (;k« 3.6) for the ECMWF computations, which suggests that the 
meridional overturning in the Atlantic thermocline is controlled by the adiabatic 
processes (advection and eddy-transfer). However, the other meteorological databases do 
not support this hypothesis. 


Database 

Adiabatic flux ratio y 
(dimensionless) 

ECMWF 

3.5675 

NCEP/NCAR 

1.1484 

SOC-GASC97 

0.3983 

Da Silva SMD94 

0.6259 


Table 6. Adiabatic flux ratios for the upper Atlantic Ocean computed from the four 
meteorological databases. 
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4. Atlantic Contribution to Global MOC 

In order to estimate the Atlantic contribution to the global thermohaline 
circulation, we can compute the diabatic and adiabatic volume fluxes over the density 
interval [24.2, 28.3] and compare them to the values corresponding to the global ocean. 

For this computation, we include in the Atlantic basin the southern ocean between 
Cape Horn and Cape of Good Hope. Note that this extension of the Atlantic surface adds 
external sources of water and might lead to uncertainties in and 1/estimates. In 

particular, we may expect differences on the order of 5-10 Sv between the entering ACC 
flow in the Drake Passage and the exiting flow to the south of Africa. Such differences 
are not significant relative to the ACC mean volume flux (100 Sv), but may be relevant 
for our results (on the order of 10 Sv). 

Figures 33 and 34 show 14(cr) and l/(cr) over the potential density anomaly 
range of [24.2, 28.3]. Note that K(cr) again represents the MOC adiabatic advection 
over the density range [cr, ], as in the global ocean computations. 

Comparisons between Figs. 33 and 25 and between Figs. 34 and 27 show many 
similarities in the sign and magnitude of the MOC components. These similarities 
suggest that the global pole-to-pole MOC is strongly influenced by the Atlantic Ocean. In 
particular, the peak values in and 17 profiles at cr« 25.35, which correspond to the 
Gulf Stream signature, are noticeable on Figs. 25, 27, 33, and 34. 

Table 7 provides the adiabatic flux ratios for the entire Atlantic basin. The very 
high value (f«4.11 for ECMWF climatology) supports the conclusion that the 
meridional overturning in the whole Atlantic Ocean is mainly adiabatic. 

Figures 32 and 34 display the discrepancies among the 1/ values obtained from 
the four climatologic databases. Again we see that the differences in 17 (cr) are not as 
large as those observed for (cr). 
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Figure 33. Diapycnal interior flux in the entire Atlantic Ocean computed from the four 
meteorological databases. 
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Figure 34. Adiabatic interior flux 1/ in the entire Atlantic Ocean computed from the four 
meteorological databases. 
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Database 

Adiabatic flux ratio y 
(dimensionless) 

ECMWE 

4.1055 

NCEP/NCAR 

2.5272 

SOC-GASC97 

1.7821 

Da Silva SMD94 

2.0053 


Table 7. Adiabatic flux ratios for the entire Atlantic Ocean computed from the four 
meteorological databases. 


5. Method Comparison 

Section III-A presents two ways for computing the residual flux at the bottom of 
the mixed-layer. The first method uses integration along the isopycnal contours [Eq. (14)] 
within the mixed-layer, while the second one computes integrals over the surface 
bounded by these contours [Eq. (21)]. 

The two methods are mathematically equivalent, but, in practice, may lead to 
relatively significant differences in estimates of 1/ and , depending on the grid 

resolution. Therefore, comparing their results gives a consistency check of the estimate 
accuracy. Note that the second method is more convenient, and for this reason it [Eq. 
(21)] will be mostly used. 

Eigure 35 gives 1^ (cr) and K (cr) profiles computed by the two methods from 
ECMWE climatology. These profiles show that and estimates weakly depend on 

the method used. This is also observed for the adiabatic flux ratios, which are, 
respectively, 7 = 3.6607 and 7 = 3.5675, for the integration along isopycnal contours 
and the surface integrals. 
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Figure 35. MOC components in the Atlantic Ocean computed from ECMWF climatology, 
(a) Diapycnal volume flux . (b) Adiabatic volume flux . 


C. PACIFIC AND INDIAN OCEANS 

The MOC components for the Indian Ocean have not been computed separately. 
Indeed, the Asian continent constitutes a northern boundary at low latitudes and therefore 
inhibits any isopycnal transport in the Indian basin, as shown by the distribution of 
potential density at 40 m depth in Fig. 9. 
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The large Pacific Ocean and its density field permit a pole-to-pole overturning 
circulation, which has to be studied. In addition, the numerous “holes” in its western 
lateral boundaries, and particularly the Indonesian archipelago, make it necessary to take 
into account the density forcing over the Indian Ocean as a possible source of water mass 
transformation in the Pacific basin. For these reasons, calculations will be carried out for 
the Pacific Ocean only and for the combined Indian-Pacific Ocean, and results compared. 

Similar to the Atlantic Ocean computations, we first limit our study of the Pacific 
and Indian Oceans to the density range for which all isopycnals are zonally bounded by 
land throughout the year: the monthly-averaged potential density distribution [not shown] 
provides the potential density anomaly range of [23.8, 25.4]. It is also interesting to 
assess the Pacific and Indian contributions to the global MOC by computing (cr) and 

(cr) over the sigma range [24.2, 28.3], used for the global ocean calculations. 

1. Diapycnal Volume Flux 

The diapycnal interior flux () in the upper Pacific Ocean is plotted in Fig. 36. 

Its average value over the potential density anomaly range [23.8, 25.4] is about 3.6 Sv 
(ECMWF climatology). This mean value is less than that computed in the upper Atlantic 
Ocean, but, at each cr, is generally higher in the Pacific than in the Atlantic. 

We note some peak values on the profile, but no characteristic signature of the 

Kuroshio Current (cr« 25.0). Indeed, we would have expected a positively high value 

around cr = 25.0 due to the high upward net heat flux at the surface of this warm western 
boundary current [see Fig. 10]. Instead, we observe a negative peak value at cr»25.1. 
This result suggests that the Gulf Stream has a much bigger impact on the MOC than the 
Kuroshio Current. In fact. Fig. 10 indicates that the highest density forcing is in the Gulf 
Stream region. 
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Figure 36. Diapycnal interior flux in the Pacific Ocean computed from the four 
meteorological databases. 


Again, we observe large discrepancies among estimates from the different 

databases, although results from ECMWF and NCEP/NCAR re-analyses show many 
similarities. Concerning SOC-GASC97 and Da Silva SMD94 climatologies, the results 
are again unrealistic, which is most likely due to the lack of observational data for the 
southern ocean. 

The profiles for the combined Indian-Pacific Ocean are close to those for the 
Pacific Ocean only. As shown in Fig. 37, the major difference is in the value of the 
diapycnal interior flux, which increases by 13 Sv (ECMWF climatology) when the Indian 
Ocean is taken into account. This implies that the Indian Ocean contribution to (cr) is 
almost constant (about 13 Sv) over the sigma range [23.8, 25.4]. 
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Figure 37. 


Comparison of the diapycnal interior flux in the Pacific Ocean and in the 


combined Indian-Pacific Ocean. The values are computed from ECMWF and 
NCEP/NCAR databases. 


2. Adiabatic Advection 

Eigure 38 shows (cr) computed from Eqs. (21), (25), and (27) using the four 
meteorological databases. In a similar way as in the upper Atlantic calculations, K 
represents here the integrated isopycnal interior flux over the potential density anomaly 

r n dV 

range!cr(,,crI. Thus, —^>0 means a northward flux, and —^<0 corresponds to a 
d(j d(j 

southward flux. 

The integrated isopycnal flux 1/, calculated from ECMWE (reference 
climatology), is negative over most of thecr range [23.8, 25.4]. This indicates that the 
residual flow in the Pacific upper thermocline is predominantly southward, with a mean 
value of -8.6 Sv. The differences in 1/ (cr) computed from the various databases are not 

as large as in (cr), which is consistent with our results for the Atlantic Ocean. 
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Contrary to the results in the Atlantic and global oceans, the comparison of the 
ispoycnal Ekman fluxes and the isopycnal residual (total) fluxes suggests that mechanical 
forcing (wind stress) is the major driving force for the MOC in the Pacific Ocean. Indeed, 
the Ekman flux dominates the residual flux, with the same sign and tendency, over almost 
the complete density range. Note that this statement is not valid for the combined Indian- 
Pacific basin. 

The computations for the combined Indian-Pacific Ocean provide profiles that 
are similar to those for the Pacific Ocean only. As shown in Pig. 39, we observe a nearly 
uniform Indian Ocean contribution to K (‘^) *^hat is about 12 Sv over the entire sigma 
range [23.8, 25.4]. 


^XlO 


—Total ECMWF 

- Total NCEP/NCAR 
Total SOC-GASC97 
Total Da Silva SMD94 

- •Ekman ECMWF 

- Ekman NCEP/NCAR 

- Ekman Da Silva SMD94 



244 24.6 24.8 

Potential density anomaly 


Figure 38. Adiabatic interior flux in the upper Pacific Ocean computed from the four 
meteorological databases. 
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Figure 39. Comparison of the adiabatic interior flux in the Pacific Ocean and in the 

combined Indian-Pacific Ocean. The values are computed from ECMWF and 
NCEP/NCAR databases. 


3. Adiabatic Flux Ratio 

Tables 8 and 9 give the adiabatic flux ratios respectively for the Pacific basin and 
the combined Indian-Pacific Ocean. This parameter is relatively low in the two oceans: 
y~2A5 (ECMWF computations) in only the Pacific Ocean and y ^\A1 (ECMWF 
computations) when we also consider the Indian basin. This suggests that the meridional 
overturning is still controlled by the adiabatic processes (advection and eddy-transfer) in 
the Pacific thermocline, whereas it is mostly diabatic in the upper Indian Ocean. 
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Database 

Adiabatic flux ratio y 
(dimensionless) 

ECMWF 

2.4542 

NCEP/NCAR 

1.4346 

SOC-GASC97 

0.4461 

Da Silva SMD94 

0.5854 


Table 8. Adiabatic flux ratios for the upper Pacific Ocean computed from the four 
meteorological databases. 


Database 

Adiabatic flux ratio / 
(dimensionless) 

ECMWE 

1.4681 

NCEP/NCAR 

0.7260 

SOC-GASC97 

0.5042 

Da Silva SMD94 

0.5683 


Table 9. Adiabatic flux ratios for the upper Indian-Pacific Ocean computed from the four 
meteorological databases. 


4. Pacific and Indian Contributions to Global MOC 

To estimate the Pacific and Indian contributions to the global thermohaline 

circulation, Figs. 40, 41, 42, and 43 display (cr) and (cr) in the Pacific and Indian- 

Pacific oceans over the density interval of [24.2, 28.3]. The chosen density range includes 
the southern ocean, where the volume flux of the ACC may result in external sources of 
water and might lead to uncertainties in and 1/estimates. 
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Figure 40. 


Diapycnal interior flux in the entire Pacific Ocean computed from the four 
meteorological databases. 
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Figure 41. Adiabatic interior flux in the entire Pacific Ocean computed from the four 
meteorological databases. 
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In Figs. 41 and 43, I4(cr) represents the integrated adiabatic advection over the 
density range [cr,cr^^]. In Fig. 41, we see that 1/ is dominated by the Ekman flux. 
Significant differences between 1/ and occur only for the low densities (cr < 25.3): 
this feature may be partially due to the external water sources in the southern ocean 
affecting the integrated flow (cr). The mechanical forcing (wind stress) also drives the 

MOC in the combined Indian-Pacific Ocean, except at low densities(cr< 25.2). The 

Ekman flux dominates the residual flux, with the same tendency, over almost all the 
density range [25.2, 28.3]. 

The comparisons between Eig. 42 and Eig. 25 and between Eig. 43 and Eig. 27 
indicate the large contribution of the density forcing over the Indian Ocean to the global 
MOC circulation (on both diabatic and adiabatic components). Indeed, the correlations 
between the global MOC components and those computed over the Pacific Ocean only 
are relatively low, whereas the same correlations with the combined Indian-Pacific Ocean 
instead of only the Pacific Ocean are much higher. This result is surprising, since the 
Indian Ocean is the smallest oceanic basin and its topography also inhibits any pole-to- 
pole MOC circulation. 

Tables 8 and 9 give the adiabatic flux ratios respectively for the whole Pacific 
basin and the whole combined Indian-Pacific Ocean. Relatively high values of f indicate 
that the adiabatic advection is also an important process in the Pacific and Indian Oceans. 
In particular, the isopycnal volume flux tends to dominate the MOC at high densities 
(southern ocean). 
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Figure 42. Comparison of the diapycnal interior flux in the entire Pacific Ocean and in 

the combined Indian-Pacific Ocean. The values are computed from ECMWF and 
NCEP/NCAR databases. 



Figure 43. Same as Fig. 42, but for the adiabatic interior flux 1/ . 
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Database 

Adiabatic flux ratio y 
(dimensionless) 

ECMWF 

2.5561 

NCEP/NCAR 

2.8388 

SOC-GASC97 

0.9889 

Da Silva SMD94 

1.6378 


Table 10. Adiabatic flux ratios for the entire Pacific Ocean computed from the four 
meteorological databases. 


Database 

Adiabatic flux ratio / 
(dimensionless) 

ECMWF 

3.0824 

NCEP/NCAR 

3.1211 

SOC-GASC97 

1.5104 

Da Silva SMD94 

2.3517 


Table 11. Adiabatic flux ratios for the entire Indian-Pacific Ocean computed from the four 
meteorological databases. 

D. CONSISTENCY OF THE ADIABATIC ADVECTION ESTIMATES 

The isopycnal flux (cr) is defined by a linear average of two net residual fluxes 

at the bottom of the mixed-layer [Eq. (25)]. Since this assumption is a critical element of 
our method for estimating the residual overturning circulation, the formulation for 1/ 
must be tested for consistency. 

One way to test 1/ estimates is to use an oceanic General Circulation Model 
(GCM). Indeed, the residual isopycnal advection diagnosed from the air-sea density 
fluxes can be compared with the values derived directly from the GCM-simulated 
velocities. 


1. Model Characteristics 

The model used for this study is based on the GFDL MOMS code (Pacanowski 


and Griffies 1999). Its horizontal resolution is two degrees in longitude and latitude. It 
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covers the whole ocean from 78°S to 84°N and is divided into 25 levels in the vertical 
direction, with its resolution increasing from 17 m at the surface layer to 510 m at the 
bottom. The bathymetry of the model is derived from the Scripps Topography. Vertical 
diffusivity varies from 0.25x10^ 'at the surface to 1.0x10 'at the bottom. 
This profile reflects the increase of the vertical mixing from the thermocline to the deep 
ocean (Bryan and Lewis 1979) and the intensification of mixing by rough bottom 
topography (Polzin et al. 1997). Heat and salt transports by mesoscale eddies are 
parameterized by the Gent-McWilliams scheme (Gent and McWilliams 1990) with 
coefficients for isopycnal diffusion of tracers and isopycnal thickness of 500 '. The 

K-profile parameterization (KPP) scheme (Large et al. 1994) is used to represent 
turbulent mixing within boundary layers. 

The simulated density field at 40 m depth is given in Fig. 44a: due to the coarse 
resolution of the model, small differences with observations remain at high latitudes of 
the North Atlantic and in the ice-covered areas. Surface density fluxes are quantitatively 
consistent with observations, and the largest differences with the observational data 
appear in regions where the coarse resolution distorts simulations of the currents (at the 
equator, in the Western Boundary Current areas, and in the Antarctic Circumpolar 
Current region). 

Heat fluxes into the ocean are determined by using conventional bulk formulas 
and are corrected for stability. Daily values for the two-meter air temperature and 
humidity, ten-meter wind speed, and zonal and meridional components of the wind stress 
have been taken from NCEP/NCAR reanalysis for the period of 1979-2001 and 
climatologic monthly means are used for all other atmospheric variables. Cloud cover 
and solar radiation are given by the International Satellite Cloud Climatology Project 
(ISCCP). Freshwater fluxes are taken from Jiang et al. (1999) and include river runoff 
data. The model is coupled to a thermodynamic sea-ice model (Visbeck et al. 1998). 

For the analysis, the 1990 simulation has been chosen: it represents the mean year 
of the studied period and does not correspond to an El Nino/Ea Nina event. The density 
flux distribution computed from the model is given in Eig. 44b. 
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The MOC estimates provided by the model are generally coherent with the 
previous studies. In addition, density fluxes, stratification, and velocity fields in GCM- 
simulation are all dynamically consistent and readily available, which enables us to verify 
main assumptions of the theory in Section 11. 

2. Adiabatic Advection 

The isopycnal flux has been estimated over the interval [24.4, 28]. This range is 
determined with the same criteria as those used in Section III-A-3: cr„ = 24.4 is the 
lightest density anomaly whose outcrop forms two non-intersecting contours of surface 
density throughout the year, and = 28 is the highest density of the ice-free sea 
surface area. 

The adiabatic MOC component 1/ (cr) has been computed in two ways: (i) by 
using Eq. (26) to infer 1/ from the sea-surface density flux, and (ii) by evaluating 1/ 

directly from the meridional velocity distribution. More precisely, for the second 
calculation, we consider a zonal section exactly at the equator, and the isopycnal 
transport is obtained by vertically integrating the zonal residual flow (the sum of the 
Eulerian and eddy-induced velocities of the Gent-McWilliams scheme) starting from 
CTmax towards lighter densities. The sign convention in both calculations is the same as in 
Section IV-A-2. 


72 




a) 



28 

27 

26 


25 


24 

23 



22 

21 


b) Surface density flux 



Figure 44. GCM-simulated fields, (a) Potential density anomaly distribution at 40 m depth. 

(b) Surface density flux from the atmosphere into the ocean (from Radko et al. 
2006). 
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For the individual basin computations, the water mass transport from other basins 
has been taken into account. In particular, the ACC contribution has been calculated as a 
divergence of the zonal residual flow within the ACC between 64°W and 18°E, which are 
the longitudes of, respectively, the Drake Passage and Africa. Over the southern Atlantic 
area, we observe a negative divergence: this suggests some inflow of water from the 
ACC into the Atlantic Ocean. This ACC convergence has been incorporated into the 
term of Eq. (25) for the Atlantic and Indian-Pacific calculations, and the corresponding 
estimate has been modified accordingly. 

Eigures 45, 46, and 47 plot 1/ (cr) for the global ocean, the Atlantic and Indian- 
Pacific basins, respectively. On these plots, we see that the 1/ profiles from the surface 
density flux (blue curves) and those obtained directly from the velocity distribution at the 
equator (green curves) match very well, especially for the global calculations. Although 
some significant differences occur at specific cr values, overall, the results in Pigs. 45, 
46, and 47 support the diagnostic model for estimating 1/ (cr), introduced in Section II-D 
[Eq. (25)]. 

Por the global ocean case, the profile of the adiabatic MOC component computed 
from observations [Pig. 27] and that based on model calculations [Pig. 45] have many 
similarities. In particular, we observe on both profiles a positive peak value of 5 Sv at 
cr« 25.7 and a negative peak value at cr« 27.1 in Pig. 27 (observations) and at cr« 27.3 
in Pig. 45 (model). However, the magnitude of (cr) computed from observations is 
generally higher than in the model. The 1/ profile given by the model does not show any 
Gulf Stream signature (negative peak value at cr« 25.35), whereas it clearly appears on 
the observational curve. This may be attributed to poor resolution of WBC Gulf 
Stream density field by the coarse-resolution GCM. 

Por the Atlantic Ocean computations, (cr) profiles from observations [Pig. 35] 

and GCM [Pig. 46] are quite similar, but, as mentioned above, they strongly differ at the 
Gulf Stream densities (around cr = 25.35). 
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K (‘^) profiles for the Indian-Pacific Ocean [Fig. 43 and Fig. 47] also show many 

similarities between the observations and the model output. Again, the magnitude is 
generally much higher for the climatologic estimates than for the GCM-based analysis. 

Concerning the wind-driven component of (cr) (Ekman isopycnal flux), the 

model outputs are characterized by a much larger contribution of this component to the 
residual (total) isopycnal flux than the observational estimates. However, both the 
observational results and the model suggest that the local wind-stress forcing does not 
drive the residual MOC adiabatic advection in the Atlantic Ocean. 


Global 



Figure 45. Adiabatic interior flux 1/ in the global ocean computed from the GCM. The 
displayed values are in sverdrups (Sv). The blue curve represents K (cr) 
computed from the simulated surface density flux. The green curve gives 
1/(cr) obtained directly from the velocity distribution at the equator. The red 

curve is the wind-driven component of K (cr) (Ekman adiabatic interior flux) 
[from Radko et al. 2006]. 


75 



Atlantic 



Figure 46. Same as Fig. 45, but for the Atlantic Ocean. 


Indo-Pacific 



Figure 47. Same as Fig. 45, but for the Indian-Pacific Ocean. 
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V. DISCUSSION AND CONCLUSIONS 


The oceanic Meridional Overturning Circulation (MOC) is accompanied by the 
diabatic transformation of water masses and interacts with the air-sea buoyancy fluxes 
(Walin 1982). The two major mechanisms for water mass transformation are the interior 
diapycnal diffusion and the diabatic processes in the upper mixed-layer. Their relative 
importance is still uncertain and much debated. The classical view (e.g., Robinson and 
Stommel 1959; Munk 1966) emphasizes the water mass transformation in the ocean 
interior and implicitly assigns it a critical role in maintenance of the meridional 
overturning. This idea was challenged by a number of modeling studies (e.g., Bryan 
1987), which demonstrate that reproducing overturning of realistic strength via diffusive 
mechanisms requires unrealistically high values of vertical mixing. Mechanical forcing 
by winds, on the other hand, appears to be an essential ingredient of the thermohaline 
system; neglecting the wind stress in multi-century simulations, for instance, leads to a 
complete shutdown of the conveyor belt circulation (Timmermann and Goosse 2004). 
These findings suggest that the interior diffusion may not be as critical for maintenance 
of the MOC as adiabatic processes (wind forcing, advection, and eddy transfer). This 
latter statement has gained considerable attention recently (Toggvieler and Samuels 1998; 
Marshall and Radko 2003; Radko 2006), but requires independent and objective support. 

The diagnostic analysis in this study may help to resolve this dispute by 
quantifying the magnitude and distribution of the isopycnal MOC component (17 ) and 

relating it to the surface forcing fields. While the isopycnal transport is not reflected in 
the interior water mass transformation, it interacts with the diabatic mixed-layer, making 
it possible to infer its pattern directly from the sea surface data. The theoretical 
framework is based on the residual mean theory that takes into account adiabatic transfer 
by the geostrophic eddies. As indicated in the schematic diagram in Fig. 6b, each density 
layer in our model can be thought of as a “leaky pipe” that is connecting the two 
hemispheres and is pumped at its two ends with different forces. It is this difference in 


77 



pumping that ultimately determines the direction and the average transport of the flow 
through the pipe, although a certain amount of water “leaks” in the interior due to 
diabatic processes. 

This study is mainly focused on the along-isopycnal volume flux , representing 
the adiabatic advection within a specific density layer. (cr) has been defined as the 

linear average of the two net residual fluxes at the bottom of the mixed-layer in the 
Northern and Southern hemispheres. Moreover, considerable attention is given to 
properties of the diapycnal flux V^^a), representing the integrated flux across the 

isopycnal surface cr in the ocean interior. Since the residual flow is assumed to be non- 
divergent, it is exactly equal to the net volume flux from the ocean interior into the 
mixed-layer within the area defined by the mixed-layer density condition cr^ < cr. 

In addition to the computations of the adiabatic and diabatic interior fluxes, the 
adiabatic flux ratio defined in Eq. (28) has been used to determine the dominant 
component of the MOC. This parameter appears to represent accurately the relative 
strength of the two components by considering the averages of the volume fluxes. 
However, it is very sensitive to the accuracy of the observational fields and also to the 
considered density range. For this reason, the values of the adiabatic flux ratio 7 must be 
interpreted with caution. 

For the Atlantic, our results show that along-isopycnal volume fluxes greatly 
exceed the corresponding diapycnal transport (7 « 4.11 for ECMWF climatology). These 
findings support the notion that the dynamics of the meridional overturning in the 
Atlantic Ocean are dominated by the adiabatic processes (advection and eddy-transfer). 
For the Pacific, the analysis suggests that the effects of the water mass transformation and 
along-isopycnal advection are mostly comparable (7 ~ 2.56 for ECMWF climatology). 
For the global ocean, the results suggest that the global thermohaline circulation is 
dominated by the adiabatic advection (7 « 3.50 for ECMWF climatology). 

In addition to quantifying the strengths of the MOC components, this study offers 
some insight into the mechanics of the pole-to-pole meridional overturning. The total 
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along-isopycnal residual circulation (1/) is divided into two distinct components: the 
mean advection also called the wind-driven component of V^, and the eddy- 

driven component. The GCM results point out that is controlled by the mean 

advection, directly associated with the Ekman pumping, everywhere except in the 
Atlantic Ocean. On the contrary, the observational results show that the Ekman advection 
is comparatively weak and not correlated with the residual adiabatic advection, except in 
the Pacific Ocean; this emphasizes the global dominance of the eddy-driven advection in 
the adiabatic interior flux. 

It is also essential to keep in mind that the application of the water mass 
transformation theory to the sea surface data is greatly restricted by the observational 
uncertainties: the air-sea fluxes in various datasets may differ by as much as a factor of 
two or more. However, this limitation is not nearly as severe for the isopycnal MOC 
component (1/) as it is for the diapycnal flux (V^). estimates are based on the 
difference between the air-sea fluxes at the Northern and Southern outcrops of each 
isopycnal surface. Therefore, for (cr) computations, our method results in a partial 

cancellation of systematic errors in flux measurements. In fact, the four climatologic 
databases (ECMWE, NCEP/NCAR, SOC-GASC97, and Da Silva SMD94) are mutually 
consistent in terms of predicting the isopycnal transport 1/ . 

Several databases have been used in this study: Eevitus 94 for the oceanic 
parameters, ECMWE, NCEP/NCAR, SOC-GASC97, and Da Silva SMD94 for the air-sea 
fluxes. This has permitted a quality control of the results and an assessment of the 
sensitivity of the calculations to the observational fields. Eurthermore, the computations 
based on GCM-simulated parameters has validated our method of estimating the 
adiabatic interior flux (V^ (cr)). 

Statistical calculations in this study offer suggestions for the database selection 
and for the optimal choice of the operational parameters in diagnostic models. In 
particular, the mixed-layer density here is represented by the potential density anomaly at 
40 m depth, and ECMWE re-analysis appears to be the best compromise regarding the 
accuracy of fluxes at the sea surface. 
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Finally, this study offers new insights into mechanics of the thermohaline 
circulation and its interaction with the atmosphere. These insights are supported by an 
agreement (Section IV) between our results and the corresponding estimates of the 
oceanic circulation from in-situ measurements and from the earlier Walin-type diagnostic 
calculations (Donners et al. 2005). 
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